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We study the half-breathing phonon in the three-band model of a high temperature superconduc- 
tor, allowing for vibrations of atoms and resulting changes of hopping parameters. Two different 
approaches are compared. From the three-band model a t-J model with phonons can be derived, 
and phonon properties can be calculated. To make contact to density functional calculations, we 
also study the three-band model in the Hartree-Fock (HF) approximation. The paramagnetic HF 
solution, appropriate for the doped cuprates, has similarities to the local-density approximation 
(LDA). However, in contrast to the LDA, the existence of an antiferromagnetic insulating solution 
for the undoped system makes it possible to study the softening of the half-breathing phonon under 
doping. We find that although the HF approximation and the t-J model give similar softenings, 
these softenings happen in quite different ways. We also find that the HF approximation gives an 
incorrect doping and q dependence for the softening and too small a width for the (half-)breathing 
phonon. 



I. INTRODUCTION 

The electron-phonon coupling for high-T c cuprates has 
recently attracted much interest. Lanzara et al^ dis- 
covered strong coupling to a mode at 70 meV in many 
cuprates. The coupling was ascribed to a half-breathing 
phonon along the (1,0,0) direction. This is an in-plane 
bond-stretching mode, where the vibrations are primar- 
ily due to two of the four O atoms surrounding a Cu 
atom in the Cu02 plane. Lanzara et alX deduced a 
rather strong apparent electron-phonon coupling A ~ 1. 
The half-breathing phonon shows an anomalous soften- 
ing when the cuprates are doped, in particular towards 
the zone boundary i 2 i 3 i 4 i 5 The softening of other phonons 
upon doping can be explained as a screening of the ions 
in the doped system^ The softening of the half-breathing 
mode, however, cannot be described in a shell model with 
conventional parameters^ This supports the idea that 
this phonon has substantial electron-phonon coupling, 
which would lead to a reduction of the frequency in the 
doped but not in the undoped system. This phonon has 
an appreciable broadening, 3 which can also be explained 
in terms of a substantial electron-phonon coupling. From 
the broadening one could estimate an electron phonon 
coupling for this mode of the order of A ~ 0.2 — 0.3, 
using the formula of Allen 6 and the density of states 
of Mattheissi Anomalous behavior of bond-stretching 
modes has also been observed in other compounds pSaLiS 

The half-breathing phonon has been studied 11-1213-14 
within the t-J model^ and a substantial softening was 
found ^Li2ii4 Phonons and the electron-phonon interac- 
tion have also been studiedi§i 1 ii& extensively within the 
local-density approximation (LDA)ii^ Bohnen et ali& 
found phonon frequencies in good agreement with ex- 
periment for YBa2Cu307. In particular, the frequency 
of the half-breathing mode along the (1,0,0) direction 
was found to be anomalously soft, in agreement with 
experiment. Since LDA does not describe antiferromag- 
netism in the undoped system, however, it cannot prop- 
erly describe phonons in the undoped case. It is therefore 



not clear how much the phonon is softened under dop- 
ing in the LDA. Furthermore, LDA calculations show 
a weak electron-phonon coupling to the half-breathing 
phonon, with A at the zone boundary being ~ 0.01. 18 
This is in disagreement with the large width of the half- 
breathing phonon, which is believed to be due to a rather 
strong electron-phonon coupling. The weak coupling 
raises questions about the reason for the low frequency 
of the half-breathing phonon in the LDA calculation for 
the doped system. 

To address these issues, we here study the three-band 
model^S including a Cu d x i_ y i and two O p-orbitals per 
unit cell. We allow displacements of the atoms from their 
equilibrium positions and take the corresponding changes 
of hopping integrals into account. In this model we can 
study phonons. We solve the model using a Hartree- 
Fock (HF) mean-field approximation, which may be ex- 
pected to simulate features of the LDA for the doped sys- 
tem. For instance, we find a similar width for the half- 
breathing phonon as in the LDA. In contrast to LDA, 
however, this approximation gives an antiferromagnetic 
solution for the undoped system. We can therefore obtain 
the softening upon doping within this framework. The 
softening of the half-breathing phonon is indeed found to 
be of the same order of magnitude as the experimental 
result, supporting the idea that the LDA can describe 
the softening. 

Alternatively, a t-J model with phonons can be derived 
from the three-band model. This model is solved using 
exact diagonalization^ Thereby, many-body effects are 
included, and the interplay between electron-phonon and 
electron-electron interactions are considered in a more 
explicit way. In this paper we compare the results in 
these two approaches. 

We find that the HF solution of the three-band model 
and the t-J model give a comparable softening of the 
half-breathing phonon. This happens, however, in a very 
different way in the two approaches. Furthermore, we 
find that the dependence on doping and on q is rather 
different in the two approaches. In the t-J model, the 
softening is ~ 6 for small 5, while the doping dependence 
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is weaker in the HF solution. The t-J model gives a 
smaller softening for the q = (1, l)n/a breathing mode 
than for the q = (1, 0)ir/a half-breathing mode, while the 
opposite is found in the HF approximation. In both cases 
t-J results are in better agreement with experiment. 

In Sec. II we introduce the models and in Sec. Ill 
the methods. The results are presented in Sec. IV. The 
results of the HF approximation for the three-band model 
are compared with the exact diagonalization results of 
the t-J model in Sec. V. 



II. MODELS 

We study the three-band model with N unit cells 
H = £d E n ida + e p n jpa + E n id ^n idi (1) 

2(7 jtr i 

+ E WfrLtp* + ^ + E 

where riida and rij pa are the occupation numbers for the 
N 3d-orbitals and the 2N 2p— orbitals, respectively, and 
ipido- and ipjpo- are the corresponding annihilation opera- 
tors. The site energies are e d and e p . There is a Coulomb 
integral U, acting between two 3d-electrons on the same 
Cu site. The model includes hopping between nearest 
neighbor Cu and O atoms and between the O atoms 
which are nearest neighbors of a particular Cu atom^i 
These hopping integrals are given by t pd and t pp , respec- 
tively, with the signs determined by the relative orienta- 
tions of the orbitals involved. Displacing an atom leads to 
changes of both electrostatic potentials and hopping inte- 
grals. Electrostatic potentials are screened differently in 
the doped and undoped systems. It was found, however, 
that the change of screening does not strongly influence 
the half-breathing phonon. 3 We therefore assume here 
that the changes of hopping integrals are the dominating 
effects^ 

We assume the hopping integral t pd has a power de- 
pendence on the atomic separation d, i.e., 

tP d {d)= (^/pj\ P m (2) 

where a is the lattice parameter. The exponent 3.5 for t pd 
was estimated from LDA band structure calculations. 23 
We neglect the phonon modulation of tpp. 

This model can be solved directly as it stands within 
some approximation, here the Hartree-Fock approxima- 
tion. We can, however, also derive a t-J model with 
phonons, and then solve this model. The t-J model has 
one site per Cu atom. In the undoped system, there is one 
hole per site. Doping introduces additional holes, which 
primarily sit on O sites. These O holes form Zhang- 
Rice singlets with neighboring Cu holes^ It is then as- 
sumed that each site is occupied by either a d-hole or 
a Zhang- Rice singlet. Expressing the atomic displace- 
ments in terms of phonon operators and working to linear 



order in the displacements leads to the t-J model with 
phononsia 

H t .j = J^( S ,.S r ^)-i^ clc ja (3) 

<i,j> <ij>c 

+E^( q ) [ b w b w + 2) 

+ E 4r g *r E 9i] J i^ + b -^)- 

ija qf 

Here c| CT creates a rf-hole on site i if this site previously 
had no hole, i.e., if it had a Zhang-Rice singlet. The 
operator 6^ v creates a phonon with wave vector q, index 
v and frequency ^(q). &i,j-x = 1 if site j is one site 
to the right of i in the x-direction. The formulas for 
<7*J (q, v) have been given elsewhere^ While we treat 
the three-band model in the electron picture, we treat 
the t-J model in the hole picture. 

In deriving this model, we assumed that t pd <^ \e d —e p \, 
U . This is a rather poor approximation for realistic pa- 
rameters. Here we have somewhat artificially reduced the 
value of t pd0 to 1.2 eV, which gives a realistio24 value of 
t = —0.47 eV if perturbation theory is used, and probably 
also more realistic values for gtj J , by compensating for 
some of the effects of perturbation theory^ The treat- 
ment above neglected quadratic terms, although these 
also can give a doping dependent contribution to the 
phonon softening. 

III. METHODS 

We use a frozen phonon method to calculate the 
phonon frequency in the HF approximation for the three- 
band model. This gives the effective Hamiltonian 

HcS = E( £ d + U (nid-a))nida + £ P E n w<y ( 4 ) 

ia jo 

-V E("^T>( n ^i) 

i 

+ E [fVLvw + h.c] + £ tgvkvw 

where {fiid-a) is the expectation value of riid-a- We also 
introduce the effective level 

ef a = e d + U{n ld - a ). (5) 

A static distortion is built into the lattice and the 
change of the total energy is calculated. From this the 
interaction with the electrons can be deduced and the 
softening of the phonon obtained. The displacements of 
the atoms in the lattice is denoted by {u}. The energy 
of the electronic system is then 

E[{u}}=Tr(H[{u}]p[{u}}), (6) 
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FIG. 1: Diagram describing the linear response, g ~ dH/du 
and g c s ~ dH c s/du describe the bare and screened pertur- 
bations due to a phonon. 



where p[{it}] is the density matrix, obtained in some 
approximation, e.g. the HF approximation. Using the 
Hellmann-Feynman theorem, it follows that 
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d 2 E 
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From the total energy calculation, we can deduce both 
d 2 E/du 2 and (d 2 E/du 2 Y 1 \ i=l, 2. {d 2 E/du 2 )^ de- 
scribes how a first order change in the external Hamilto- 
nian leads to a first order change in the density matrix, 
which acts back at the Hamiltonian. This contribution 
can be calculated in linear response 



dp 
dv 



x- 



dH c 



dv 



(8) 



where \ is the response function for nonintcracting elec- 
trons and the derivative of the effective Hamiltonian is 
obtained from the total energy calculation. Thus we have 



d 2 E 
du 2 



(i) 



= Tr 



— ] = TV 

du du 



dH_ ggeff 

du du 



We introduce a dielectric function 



Tr (—X—) /Tr ( dH x dH <* 
\du du J \du du 



(9) 



(10) 



The result in Eq. © can also be obtained from a di- 
agrammatic technique. It is described by a bubble di- 
agram, as shown in Fig. ^ This is the only diagram 
which enters in the HF approximation. We observe that 
the screening of the perturbation, described by g c g, en- 
ters at one of the vertices. Using g c g at both vertices 
would lead to double counting. By using the diagram- 
matic approach, we can obtain not only a contribution 
to the energy but also the width of the phonon. 



IV. RESULTS 

Below we present results for the three-band model in 

the Hartree-Fock approximation, using the parameters 
t pdo = L6 eVj t p P = o and jj = 8 eV _ The i attice pa- 
rameter is a — 3.8 A. We have adjusted e p so that the 
separation between the effective 3<i-level and 2p-levels 
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FIG. 2: Softening Au) p h for the zone boundary half-breathing 
phonon. Results are given in the HF approximation (full 
line), the t-J model (dashed line) and according to experiment 
(dash-dotted line) as a function of the hole doping S. The lines 
serve as a guide for the eye between the few points in the t-J 
model and according to experiment. The HF approximation 
refers to results for the shift in a paramagnetic calculation 
for doping S relative to an antiferromagnetic calculation for 
5 — 0. The dashed part of the HF line indicates schematically 
that the systems goes antiferromagnetic for small dopings. 
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3 eV, 



(11) 



where {rider) is the average occupation of the 3d-levels 
per spin. This is a typical LDA estimate for e d s — s p s in 
a three-band model»2i We perform two calculations, one 
for the undistorted lattice and one for a lattice where a 
phonon has been built in. This gives the second deriva- 
tive, d 2 E/du 2 . From this derivative we obtain the soft- 
ening of the phonon due to the interaction with the elec- 
trons in the model, reducing the frequency u> p ho to to p h- 
Here u> p ho is assumed to be due to forces not included in 
the model in Eq. Q, e.g., electrostatic forces and core- 
core overlap effects. The calculations were performed for 
a cluster of 32 x 32 CuC>2 units and periodic boundary 
conditions. The doping was chosen in such a way that de- 
generate levels were either completely full or completely 
empty, i.e., all "shells" were either full or empty. 

We first consider the half-breathing phonon for q = 
(1, 0)7r/a, where the two O atoms at (a/2,0) and (-a/2,0) 
move towards (or away from) from the Cu atom at 
(0,0). We perform a calculation for the undoped sys- 
tem, having five electrons per unit cell, and allowing for 
spin-polarization. We adjust Lu p ho so that the softened 
frequency u> p h is 0.080 eV for the zone boundary half- 
breathing phonon, as is found experimentally. 

The result as a function of the hole doping 6 (6 < 
means electron doping) is shown in Fig. [21 The dop- 
ing dependence is relatively weak. Since the undoped 
system is antiferromagnetic but the doped system is (as- 
sumed to be) paramagnetic, as found experimentally, the 
softening is not necessarily small for small dopings. The 
spin-polarized system has a large gap of about 4.6 eV. 
Due to this gap, the response of the system to a phonon 
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is substantially weaker than for the paramagnetic state. 
For instance, [d 2 E/du 2 ]^ contributes a softening of only 
about 4 meV for the spin-polarized system but 6 meV 
for a non spin-polarized system with a similar number 
of electrons. Similarly, the contributions of [d 2 E / du 2 ]^ 
are about 8 and 13 meV, respectively. 

We find that screening reduces the quantity 
[d 2 E/du 2 ]^ by about a factor of two, i.e., e w 2 
in Eq. i|10|) . By evaluating the diagram in Fig. ^ we 
find that the zone boundary half-breathing phonon is 
broadened by about 0.4 meV for 6 = 0.16. This is similar 
to what Bohnen et alJ& found in an LDA calculation. 

It is interesting to compare the half-breathing phonon 
and the q = (1, l)ir/a breathing phonon. In the latter 
case all four O atoms surrounding a Cu atom move to- 
wards this atom. In the model considered above (t pp = 0) 
and for 6 = 0, the model has nesting for q = (1, l)-7r/a. 
The result would be a very strong response for small 5, 
making the calculation rather meaningless. Instead we 
consider t pp = 1.1 eViSi This changes the Fermi sur- 
face and avoids an unrealisticly strong nesting. In this 
case we adjust u p ho so that uj p h = 0.090 eV in the spin- 
polarized calculation for the undoped system, as found 
experimentally. For the doping 5 = 0.096 we find that 
the q = (l,l)7r/a breathing phonon is softened by 11 
meV. This is a larger softening than was found in Fig. [21 
for this doping. We have also performed a calculation for 
the q = (l,0)7r/a half-breathing phonon in this model 
and for this doping. The softening is 7 meV, which again 
is smaller than for the the breathing phonon. This is 
contrary to experiment, where a larger softening is found 
for the half-breathing phonon. 



V. COMPARISON OF THE THE t-J MODEL 
AND HF SOLUTION OF THE THREE-BAND 
MODEL 



To make contact with the t-J model, we transform 
the three-band model to a one-band model. We assume 
that the O 2p-levels are far below the Cu 3d-levels (in 
the electron picture). We can then project out the O 2p- 
levels and obtain a model with just effective Cu 3d-levels. 
This model can be compared with the t-J model, since 
both models have one band. Here we focus on the linear 
response term (d 2 E/du 2 )^ in Eq. 10, which are of the 
same order of magnitude in the t-J model and the HF 
approximation for a typical 5 ~ 0.1. 

We introduce a projection operator P, which projects 
out the Cu 3<i-levels, and its complement Q = 1 — P. 
Considering the resolvent operator of the HF Hamilto- 
nian, we obtain 



P{z-H)- 1 P 

= [P(z - H)P - PHQ(z - QHQ)~ 1 QPH]~ 



(12) 



where z is some typical energy. We then obtain the ef- 



fective one-band parameters 
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where the sum over j for e^" c runs over the nearest O 
neighbors of Cu and j for i° no refers to the common 

nearest neighbor O atom of the Cu atoms i and i . We 
first consider the unperturbed (no phonon) system. We 
choose z to be in the middle of the band and solve the 
self-consistent equations 
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(14) 
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where the factor 4 for s^ ne comes from the four O neigh- 
bors of a Cu atom. 

We then introduce a phonon in the system and ask for 
the linear response of the electronic system. A perturba- 
tion term is introduced in the one-band Hamiltonian 



,pd 



One 



(15) 



where dt pd are the changes in the Cu-0 hopping integrals. 
To linear order there is no change in i°, no . 

For a half-breathing phonon at the zone boundary, the 
on-site perturbation is 



OH 

du 



Mr, 



dt 



pd 



dr £° ne 



-u. 



(16) 



where u is the absolute value of the phonon ampli- 
tude. The quantity dH/du enters in the calculation of 
(d 2 E/du 2 )W in Eq. ©. 
In the t-J model, the on-site perturbation is 14 
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u, 



(17) 



where A = 0.96. The first term comes from the hopping of 
a 3<i-hole into the O 2p-states and the second term from 
the hopping of a O 2p-hole into the Cu 3<i-state. The 
second term has no correspondence in Eq. (|16|l . Equa- 
tion (|17l) has an additional factor 2 coming from a phase 
coherence factor in the Zhang-Rice singlet. This results 
from the singlet being explicitely written as a sum of two 
terms. Both these effects are genuine many-body effects. 
The —1 in the first term in Eq. (|17f) results from taking 
the difference in the energy gain of a Zhang-Rice singlet 
and a single 3c?-hole. 
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TABLE I: Contributions to the phonon softening in the t-J 
model and in the HF solution of the three-band model. Here 
A : 



Source 


t-J 


HF 


Ratio 


Coupling 


[(2A 2 - 1)/A + 2X 2 /(U - 


-A)] 2 (1/A) 2 


ps 3 


Sum rule 


w 2SttN 


ps tyN 


ps 25 


Screening 


1 


ps 0.5 


ps 2 


Denominator 






PS 1 


Product 






125 



We find that this on-site coupling dominates over the 
off-site coupling. 14 In the qualitative discussion we there- 
fore only consider the on-site coupling gll J (q) = ff(q). 
The phonon self-energy in the t-J model is 

s(q) 2 x*" J (q>w) 



n(q,w) = 



l + 5(q) 2 X*- J (q,w)D (q J w)' 



(18) 



where Do(q,u/) is the noninteracting phonon Green's 
function. We hnd that the second term in the denom- 
inator is small for the parameters considered here. In 
the formal discussions below we therefore neglect it. The 
phonon self-energy is then proportional to the response 
function. This result can then be directly compared with 
the result in Eq. for the one-band model. In both 
cases the response function is multiplied by the appro- 
priate coupling constant squared, given by Eq. (|16J) and 
Eq. (TJJ, respectively. 

To compare the t-J model with the one-band model, 
we put |e£ nc -£ p | = \e d -e p \ = A = 3 eV and U = 8 cV. 25 
We find that the square of Eq. (|17|) is about a factor of 
three larger than the square of Eq. 1)16(1. This difference 
due to differences in coupling constants is shown in Table 

ID 

The linear response of the one-band model is given by 



,HF 



X 



(q,c) = 2^[l-/(k + q)]/(k) 



(19) 



uj - uj(k, q) + i0+ 



■w(k,q) 



where /(k) is the Fermi function for a state with wave 
vector k and energy e(k) and w(k, q) = e(k + q) — e(k). 
We have the sum-rule 

1 f°° 

— ^2 | Im X HF (q, ^ + i0 + ) | = 4tto(1 - n)N « ttN, 

N 

(20) 

where N is the number of sites and n is the fractional 
filling of the band. Typically, we are interested in a sys- 
tem with n = (1 — S)/2 ps 0.5 electrons per site and spin, 
which leads to the right hand side of Eq. (|20[) . For the 
t-J model we ask for the carrier-carrier response function 



x'~ J (qw) 



E 



KHp(q)|o)| 5 



(21) 
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UJ — U) v + t0+ UJ + UJ V — j0+ 



where \v) is an excited many-electron state with the ex- 
citation energy uj v and p(q) is carrier density operator. 
One finds 1 ^ 



-y 

N ^ 



q#0 



Imx'" J (q,w + iO + ) = 2wN6(l - S). (22) 



Eqs. (|20|) and (|22|) differ by approximately a factor of 25 
for small 5 (see Table QJ. 

The two approaches further differ by the screening in 
the HF approach, discussed in Eq. (|10fl . This reduces 
the HF result by roughly a factor of two. Finally, we 
have to consider that a typical energy denominator in 
Eq. PU)| and Eq. J2U are different. Calculating the 
average (l/u>), we find comparable results for the two 
models. As can be seen from Table H] the difference be- 
tween the two approaches results in a ratio of about 12<5. 
For 6 ~ 0.1, the two results are then similar. This in- 
dicates why the t-J model and the HF solution of the 
three-band model can give similar softening of the half- 
breathing mode although the physics is quite different. 

We next consider the imaginary part of the phonon 
self-energy n, which gives the phonon width. As an 
orientation, we first consider a simple model. Since Im 
x(w) ~ to for small lu, we assume 



Im n(q, uj) 



Auj, 
0, 



if M < W; 
otherwise, 



(23) 



where A is some constant. From the Kramers-Kronig 
relation, we can then derive 
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pit 



(24) 



where 7 = 21m n(q, u>) is the full width at half maximum 
of the phonon and Auj p h is its shift. 

Figure [31 compares Im n(k, oj) for the half-breathing 
phonon in the one- and three-band models and the 
breathing phonon in the three-band model. The one- 
band model was constructed to describe what happens 
close to Ep and therefore the one- and three-band mod- 
els agree very well for small to. For 3> 0, the one-band 
model gives a larger |Im II(q, u>)\ than the three-band 
model, and |Re n(q, uj) \ is overestimated correspondingly 
in the one-band model. Appropriate numbers for the 
the half-breathing phonon and the one-band model are 
W = 2 eV, ujph = 0.07 eV, Auj ph = -0.007 eV. Insert- 
ing this in Eq. I|24|l leads to 7 = 0.8 meV. This is about 
twice the width actually calculated. The reason for this 
overestimate is that Im n(q, u>) is actually smaller than 
assumed in Eq. I|22l) for small uj. 

Figure[3]also shows the HF result for the q = (1, l)7r/a 
breathing phonon. These results were obtained for t pp — 
1.1 eV. There is then no strong nesting of the Fermi sur- 
face. The HF approximation, nevertheless, gives a larger 
broadening for the q = (1, l)7r/a breathing phonon than 
for the q = (1, 0)7r/a half-breathing phonon, whether the 
latter is calculated for = or = 1.1 eV. As is 
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FIG. 3: Im II(q, uf) for the q = (1, 0)n/a half-breathing mode 
in the three-band model (full line) and the one-band model 
(dashed line) and for the q = (1, l)n/a breathing mode (dot- 
ted line). The doping is 8 — 0.1. 



well-known,- 15 this is due to the fact that the wave vector 
q = (1, l)ir/a fits better to the traces of nesting left over 
for t pp = 1.1 eV. 

To study the phonon width in the t-J model, we cal- 
culate the phonon self-energy^ We calculate the phonon 
spectral function £>(q, oj) using exact diagonalization. 
A Hilbert transform can then be used to obtain the 
phonon Green's function D(q, lu). The phonon self- 
energy Il(q, lu) is calculated by inverting 



Z?- 1 (q,^) = J D - 1 (q,^)-n(q,a;), 



(25) 



where Do(q,lu) is the noninteracting phonon Green's 
function. This approach has important advantages for 
small systems. B(q,cu) has too few structure to deter- 
mine the phonon width. A broadened version of II(q, lu), 
however, can give such information^ Results for II(q, lu) 
are shown in Fig. 

In the view of Eq. (|24|l . one might have expected the 
width of the half-breathing phonon to be similar in the 
HF approximation and the t-J model, since the shifts are 
similar. This is not true, however, since the frequency 
dependence differs strongly from the linear dependence 
assumed in Eq. This is illustrated in Fig. ^ The 

figure shows that for the t-J model with J = 0, some 
spectral weight has been moved to small frequencies due 
to the hopping constraint, which creates low energy ex- 
citations. This is even more true for the finite J case 
(J/t = 0.3). The present clusters are too small to give 
reliable results, in particular for the low-lying excitations, 
and the results above should be considered as qualitative. 
They illustrate the general trend, however, of transferring 



some spectral weight to small frequencies, which tends 
to lead to a substantially increased width of the half- 
breathing phonon. 
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FIG. 4: Im II(q, to) for the q = (1, 0)n/a half-breathing mode 
in the t-J model for J/t = 0.3 (full line) and J/t — (dashed 
line) and in the one-band model (dotted line) . The results for 
the t-J model were obtained for a 4 x 4 cluster. The doping 
is 5 = 0.125. 



VI. SUMMARY 



We have studied the properties of the half-breathing 
phonon in the three-band model of a high-T c cuprate. 
The results in a t-J model with phonons and the HF 
approximation of the three-band model were compared. 
Although the two approaches give similar softenings at 
typical dopings 6, the underlying physics is quite differ- 
ent. The hopping constraint in the t-J model, resulting 
from the strong Coulomb repulsion, leads to a strong re- 
duction of the response to a phonon. This reduction is, 
however, partly compensated by several other effects. In 
particular, the coherent hopping in the formation of the 
Zhang- Rice singlet tends to increase the coupling in the t- 
J model. We find that the doping and q dependences for 
the (half-)breathing mode are quite different in the two 
approaches, with the t-J model giving better agreement 
with experiment. The HF approximation also gives a too 
small width for the half-breathing phonon. We therefore 
conclude that many-body effects play an important role 
for the interactions of the half-breathing phonon with the 
electrons. 
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P. Horsch, O. Jepsen, N. Nagaosa and Z.-X. Shen for 
many useful discussions. 
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